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ratio is 8 per cent, at ordinary temperatures which remains constant till 
about 120°, and then falls off, rapidly at first and then slowly to about 
1*2 per cent, at the critical point. There is no change of imperfection of 
polarisation on passing through the critical point. The correction due to this 
in the expression for the intensity of the scattered light is given. 

In conclusion, I have great pleasure in recording my indebtedness to 
Prof. C. V. Baman who suggested the research, and who continued to take an 
inspiring interest in its progress. The experimental work was carried out in 
the Physical Laboratory of the Indian Association for the Cultivation of 
Science, Calcutta. 
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S 1. Introduction. 

In both physical and biological science, we are often concerned with the 
properties of a fluid, or plasma, in which small particles or corpuscles are 
suspended and carried about by the motion of the fluid. The presence of the 
particles will influence the properties of the suspension in bulk, and, in 
particular, its viscosity will be increased. The most complete mathematical 
treatment of the problem, from this point of view, has been that given by 
Einstein,* who considered the case of spherical particles and gave a simple 
formula for the increase in the viscosity. We have extended this work to 
the case of particles of ellipsoidal shape. 

The second section of the paper is occupied with the requisite solution of 
the equations of motion of the fluid. The problem of the motion of a viscous 
fluid, due to an ellipsoid moving through it with a small velocity of transla- 
tion in a direction parallel to one of its axes, has been solved by Oberbeck,f 
and the corresponding problem for an ellipsoid rotating about one of its axes 
by Edwards.J In both cases the equations of motion are approximated by 
neglecting the terms involving the squares of the velocities. It may be seen, 

* " Eine neue Bestimmung der MolekUldimensionen, 5 ' ' Ann. d. Physik,' vol. 19, p. 289 
(1896) ; with a correction in vol. 34, p. 591 (1911). 
t * Grelle,' vol. 81 (1876). 
J 'Quart. Jour. Math., 5 vol. 26 (1892). 
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a posteriori, that the condition for the validity of this approximation is that 
the product of the velocity of the ellipsoid by its linear dimensions shall be 
small compared with the " kinematic coefficient of viscosity " of the fluid, 
In relation to our present problem, it will therefore be satisfied either for 
sufficiently slow motions, or for sufficiently small particles. 

The third section is devoted to the consideration of the forces acting upon 
the particle as a whole. It is found that these reduce to two couples, one 
tending to make the particle adopt the same rotation as the surrounding fluid, 
and the other tending to set the particle with its axes parallel to the 
principal axes of distortion of the surrounding fluid. 

In the fourth section we investigate the motion of a particle which is 
subject to no forces except those arising from the pressure of: the fluid on its 
surface. It is found that there is a variety of possible motions corresponding 
to different initial conditions, They are periodic and of a general character 
somewhat like the motion of a top. There is thus a certain degree of inde- 
terminacy in the problem as treated by the approximation which neglects the 
squares of the velocities. 

Sections five and six are concerned with the dissipation of energy and with 
the viscosity of suspensions. It is found that the increase of viscosity may 
still be represented by a formula of Einstein's type with a modified numerical 
factor. Owing to the indeterminacy revealed in section four, it is not possible 
to give a definite value for this numerical factor. It is, however, possible to 
specify an upper and a lower limit for this factor, and the difference between 
these limits diminishes to zero as the particles approach a spherical shape. 

Section seven is devoted to a tentative suggestion for the removal of this 
indeterminacy. It is, that when a variety of motions is possible for given 
boundary conditions under the approximated equations of motion, the actual 
motion will be that which corresponds to minimum dissipation of energy. If 
this hypothesis be granted, it is possible to state many of the results of the 
paper with greater precision. 

§ 2. The Solution of the Equations of Motion , 

Consider, in the first instance, the modification produced in the motion of 
a viscous fluid by the presence of a single ellipsoidal particle. We shall 
assume that,, apart from the disturbance produced in the immediate neighbour- 
hood of the particle, the motion is steady, and varies in space on a scale which 
is large compared with the dimensions of the particle. It appears from the 
work of Oberbeck that, if the fluid at a distance from the particle has a 
velocity of translation relative to the particle, the latter is acted upon by a 
resultant force. It will appear in the course of the investigation that any 
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other motion of the fluid can only give rise to a resultant couple. The 
particle will therefore ultimately assume the velocity of translation appropriate 
to that part of the fluid which it displaces. On the above assumptions this 
velocity of translation will be sensibly uniform, and we may therefore take 
the centre of the particle as at rest. 

Let x, y, z be rectangular co-ordinates referred to axes fixed in the particle 
and moving with it, and let the surface of the particle be the ellipsoid 

x 2 /a 2 + y 2 Jb 2 + z 2 /c 2 = 1. (1) 

The undisturbed motion of the fluid in the neighbourhood of the particle is 
then given by 

M = ax + lly + gz + rjZ—ty 

vq = hx + by + fc-f gfc- %y V, (2) 

Wo = gx + fy+Cz+Zy — rjx ^ 

where a, b, C, f, g, h, £, i], fare the components of distortion and rotation of 
the fluid. These are taken to be constant in space through a volume which 
is large compared with the dimensions of the particle. They will, however, 
vary with the time, since the particle will rotate under the influence of the 
fluid, and the axes of x, y, z will be rotating axes. 

Let the spins of the ellipsoid about its axes be <o h 6> 2 , g>3. Take a set of 
axes, x', y', z\ in fixed directions with their origin at the centre of the particle, 
and let the components of distortion and rotation of the fluid in its 

undisturbed motion referred to these axes be a', b', c', These are 

constant, but the corresponding unaccented quantities are linear functions of 
these, with coefficients which depend upon the direction cosines of the two 
sets of axes, and which therefore vary with the time in consequence of 

&>1, (02, «3- 

Neglecting squares and products of the velocities, the equations of motion 
of an incompressible viscous fluid referred to moving axes are of the type 

,u-V%— -^ = pJ _ a^ -|_ W2W , I (3) 



ix 



\dt 



where p, fi,p are respectively the density, coefficient of viscosity, and mean 
pressure of the fluid. 

The spins o>i, w 2j «3 are produced by the motion of the fluid, and, as will 
appear explicitly later, they are of the same order as the velocities. Hence to 
our order of approximation all products like (03V may be neglected. It is not 
quite so obvious that the remaining terms on the right-hand side of (3) are of 
the second order. It will appear later that u, v, w are homogeneous linear 
functions of a, b, C, ..., with constant coefficients. They are, therefore, 
homogeneous linear functions of the constants a', b', c', ..., with coefficients 
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which are functions of the direction cosines of the moving axes, Hence the 
only terms which can arise in du/dt are of the type o^a', which is of the 
second order of small quantities, and is therefore to be neglected. 
The equations of motion therefore reduce to 



ox 



O'Y) 
fl^ 2 V = 7 " 



a?/' 



jjbSjhv 



3s' 



(4: 



with the equation of continuity 



3w , 3y , 3z/; A 



3# 3// a 

We have to find a solution of (4) and (5) which agrees with (2) at great 
distances from the origin, and which gives on the surface of the ellipsoid 



jt 



u = 0)2 z — a)3y, v = co&x — coiz, 
tls usual, let X denote the positive root of 
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(D\y — co%%. 
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and 
and let 



a 2 + A. 6 2 + X c a + X 
A = {(a 2 + X)(b 2 + X)(c 2 + \)}K 



(7) 
8) 
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(a 2 -f-A)A' 
We will also write 



J3 = 



•jo 



rf\ 



L V (& 2 + X)A ; 
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x (c 2 + \)A 
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\d\ 



K (b 2 + \)(c 2 + X)A' "* _ J A (6 2 + X)(c 2 + X)A 
with symmetrical integrals for /3', /3", 7', 7", so that 



, 7-/8 
b 2 —c 
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h 2 — c' 



(11) 



We will denote the corresponding integrals, in which the lower limit of 
integration has been replaced by zero, by a , #>'> — 

We have the following well-known solutions of Laplace's equation :• — 



O = 



'00 



A. 



X 



,2 






?/ 



+ 






1 "^ d\ 



a 3 +.\" t ".& a + \ r c 2 + X J A 



12) 



and 



%i = *yz, %2 = /3'sas, %3 = 7^7/- (13> 

Ve shall seek a solution in which u, v, w are expressed as linear functions 
of the first and second differential coefficients of O and the first differential 
coefficients of %i, %2, %3- These latter can, of course, be expressed in terms of 
the second differential coefficients of O, but as the corresponding terms enter 
into the solution in a somewhat different way, it will be convenient to keep 
them separate. 
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From (7) we have 

3\,_ 2a;P 2 dx ^ 2yP 2 dx _ 2zP 2 

dx a*+X' dy b* + X' dz c 2 + \ 
where 



(14) 



?/ z 



P 2 (a 2 +'\) a (6 2 + \) 2 (c 2 ~f X) 2 
We then have a number of formulae of which the following are typical : 

and 

dx dy dz A 3 ' 

dy " {b* + Xf (c* + X) A' 9s ~ " (6» + \)(c» + X)»A' V ' 

All these differential coefficients tend to zero at infinity. 
We assume 

u = « + |(R Xl + S Xa + T X 8) + W^-V^ 

| A / 3 2 n 8n\ | H / dm dn\ Ql i dm an 

V 3-c 2 3as/ \ dxdy dy J \ 3»3« 3s 



\ 3# 2 3# 3 # dx dz} 

\ 3-^c 2 3^3?/ ' dxdzj 
v = v, + |-(E %1 + S %2 + T X3 ) + U^i-W^ 

\ ox ay oif oyoz/ 

+ Vdxhi dx) + V By* dy) + Vdydz dz I 



(18) 



and 



3^ 3y 3?/ 2 3^ dzj 
^o + |(Exx + S X2 + T % 3) + V^-U^ 



(19) 



\ d^c d3 d^ ds d* 



,2 



[ a , 3 2 n 3n\ | F / dm 3fl\ ( c ^ dm do,) (2Q) 

dxdz dx) \ 3y9z 3#/ V. 3s 3 dz f 
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may be seen by inspection that these values of u 9 v, w satisfy (5), since 
a + b + C = for an incompressible fluid, and O and %i, % 2 , %3 satisfy 
Laplace's equation. 

Substituting these expressions into (4), it is easily seen that these are 
satisfied provided that 



i ,= 1 , + 2 A t|AS + BS + o2§-f(F + F / )^ 



L ox* 



r 



dz 2 



dyd. 



] z 



+ ^ +G '>H^ H+H '>S}' < 21 > 



dx dy 

where p is the constant mean pressure at a distance from the ellipsoid. It 
therefore remains to satisfy the conditions (6) on the surface of the ellipsoid 
X = 0. Substituting from (16) and (17) into (18), (19), (20), we have 

-hy{h-Z+yT-2f3II + 2otW} 



2dP 2 



(a 2 -f- A,) A 



{E + 2 (b 2 -f X) F -f 2 (c 2 4- X) F f }yz/(b 2 + X) (c 2 + X) 



+ {S + 2(c 2 H-X)G + 2(a 2 + X)G / }^/(c 2 + X)(a 2 -fX) 
+ {T + 2(a 2 + X)H + 2(& 2 + X)H ; }^/(a 2 + X)(& 2 + X) 
+ {W-~2(a 2 -fX)A-f2(& 2 + X)B}2/ 2 /(^ 2 + X) 2 



{V-~2(c 2 + X)C + 2(a 2 + X)A}^ 2 /(c 2 H-X) i 



= %{h + Z+y'T + 2/3Il-2uW} 

+ y{b-r a / U~7 / W~2(«4-/3-|-7)B} 

4 

+ «{f_£ +a 'R_2 7 F + 2/SF} 



2</P 2 



+ X)A 



{ K + 2 (6 2 + X) F + 2 (c 2 + X) F } ^s/(& 2 + X) (c 2 + X) 



+ { S + 2 (c 2 + X) G 4- 2 (a 2 + X) G' } kc/(c 2 + X) (a 2 + X) 
4{T + 2(a 2 + X)H + 2{6 2 + X)H'}^/(a 2 + X)(& 2 + X) 
+ {U-2 (6 2 + X) B + 2 (c 2 + X) C}z 2 /(c 2 + X) 2 



(22) 



-{W-2(a 2 + X)A + 2(F + X)B}ic 2 /(a 2 + X) s 



(23) 
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and 

w — «{g— i? + £'S— 2«G+2yG'} 

+ 2/{f+£+a'B+2yF — 2/3F'} 
+ s{c + /3'V— «'U— 2(a+/3 + y)C} 
2zP 2 



(c 2 + X) A 



{R + 2(& 2 + X)F+2(c 2 + X)F'}^/(& 2 + X)(e 2 + X) 

+ {S + 2(c 2 + X)G+2(a 2 + X)G'}^/(c 2 +X)(« 2 + X) 
+ {T + 2(a 2 +X)H + 2(& 2 + X)H'}a?2//(a 2 + X)(J 2 + X) 

+ {V-2(fl 2 + X)C + 2(a 2 + X)A}« 2 /(a 2 + X) 2 

-{U-2(& 2 + X)B + 2(c 2 + X)C}^ 2 /(S 2 + X) 2 



(24) 



Putting X = 0, identifying with (6), and equating coefficients, we obtain 
equations which determine the constants. On solution we find their values 
to be as follows : — 

A = i{2ao ,/ a-iSo // b-7o ,/ c}/(/3o / V / + 7o / V / + «o /, /8o /, )' 
B = ^{2/3o /, b~ 7 o ,/ c-«o ,/ a}/(^v , +7o / V , + «o / W / ) 

C = i{2 7 o ,, C~ao ,, a-/9o ,, b}/(A , V , + 7o // «o ,, + «o // /3o ,/ ). 



>> (25) 



■p = /3 f--e 2 fto(gT-ft>]) 
2«o 7 (^o + ^7o)' 

q. = 7og~a 2 /3o(??~-ft)2) 
2/3 ' (c 2 7o + a 2 a ) ? 

■o- «oh--5 2 7o(g > — o) 3 ) 
27 '(« a *o + & 2 /9o)' 



■p, = Tof-h h 2 *o(l;-- coi) 
2a , (& 2 /3 + c 2 7o) 

_ «og 4- c 2 /3 (17 — o) 2 ) 



~\ 



G 



2/3o' (c 2 7o + a 2 «o) 

it/ _ /3oh + <x 2 7o(£— ws) 
27o'(« 2 «o + S 2 A) 



S 



(26) 



J 



E = -f/«o', S = -g//? ', T = -h/ 7o '. (27) 

The values of U, V, W will not be needed for our purpose. They are 

given by 

U = 26 2 B-2c 2 (T 

V = 2c 2 C-2a 2 A >. (28) 

W = 2a 2 A-2& 2 B_ 

Thus the constants in (18), (19), and (20) are all uniquely determined and 
the velocity of the fluid is known at all points. 
For later use we note the following relations : — 

A + B + C = 0, (29) 

U + V+W = 0, (30) 

and, using the relation a + b + C = 0, 

Aa + Bb + Cc = J («o"a 2 + /3o"b 2 + 7o"C 2 )/( y e 'V' + 7o"«o" + *o'W)- (31) 
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§ 3. The Resultant Couple on the Particle, 

For an incompressible fluid the stresses are given, with the usual nota- 
tion, by 

xx = — ^ + 2/x 7T-, y# = ^ + _ . (o J) 

The calculation of the values of these stresses on the surface of the 
ellipsoid is greatly simplified by noting that most of the brackets in (22), 
(23), and (24) vanish for X = 0, and that, therefore, the only terms surviving 
in the surface values of the differential coefficients of the velocity compo- 
nents are those arising from the differentiation of the brackets. 

If we denote the components of force per unit area acting on the surface of 
the ellipsoid by X, Y, Z, we have 

X = {xx x/a? -f xy y/b 2 + xz z/c 2 } 

Y = {yx xja 2 -f yy y/b 2 -f yz z/c 2 } Y . (33) 

Z — { zx xja 2 ~f zy y/b 2 -f zz z/c 2 } 

The reduction is straightforward and, after simplifying as far as possible by 
means of the relations between the constants, we obtain 

Y = ~^ l + ^{ H ' 5 + B l + F i} ~ W«oA + /3oB + yo C) I 

Z = _^p*+§g|G J+F' |+o||-4/*P(«,,A + /3oB+7oC)|. 

(34) 

Since, apart from P, which is even in x, y, z, only first powers of the 

co-ordinates appear in these expressions, they will separately vanish on 

integration over the surface of the ellipsoid. Hence there is no resultant 

force acting on the particle. Denoting the components of the resultant 

couple acting on the particle by L, M, 1ST, we have on integration over the 

surface of the ellipsoid 

L = JJ(jrZ-«Y)dS, 

which on substitution from (34), omitting terms which obviously vanish on 
integration, gives 



Lzr .^f r/ *» » 2X 



abc I 



Aj 
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This integral is readily evaluated in terms of the volume of the ellipsoid, 
and we have from this and two similar expressions 

7T/x (G' -G), N = Sf. itfi (H' -H). (35) 



3JL 
3 



L = ~ 3 3 2 -7r/*(F'-F), M = 
Substituting from (26), these give 



L = 



M = 



__ 167T/JL 



3(J 2 /3 + c 2 7o) 



{(6»-C»)f+ (J 8 +<!»)(? -6)0}.. 



~N 



167T/L6 



3 (c 2 y -f a 2 u ) 



{(c'-a^g + ^ + a 2 )^-^)}, > 



-pr _ 16tt/x 



^^_ { («.-J»)h + (a« + ^)(t-« 8 )}.J 



(36) 



The resultant couple acting upon the ellipsoid may thus be regarded as 
consisting of two parts : one which vanishes when the ellipsoid has the same 
resultant spin as the fluid, and one which vanishes when the axes of the 
ellipsoid coincide with the principal axes of distortion of the fluid. 



§ 4. The Motion of the Particle — Laminar Motion. 

If the particle is subject to no forces except those exerted by the fluid 
upon its surface, then in the slow motions to which our assumptions have 
already restricted us, the resultant couple on the particle must vanish at 
every instant. Hence, from (36), 

(V + c^co! = ^(£-f-f) + c 2 (£-f) - 

(c 2 + a 2 ) a) 2 = c 2 ( V + g) -j- a 2 (vS) k (37) 

(a 2 + b 2 )co s = a 2 (£ + h) + & 2 (f-h)^ 
or, from (35), 

F' = F, G' = G, H' = H. (38) 

These relations will simplify many of our formulae and, in particular, we 
may note that (26) are replaced by 

f „ * g 



F = 



2*o , (& 2 -M 2 )' 



G = 



2/% / (c 2 + a 2 )' 



H = 



(39 



2 7 o> 2 -h& 2 )* 

It has not been found possible to obtain a solution of equations (37) which 
would give the motion of the particle for the most general possible undis- 
turbed motion of the fluid. Some light may, however, be thrown on the 
problem by the consideration of a particular, but important, case. 

Take axes x', y f , z\ fixed in direction, and let the components of the fluid 

velocity in the undisturbed motion be uq\ Vq, iv ' } parallel to these axes. 

Consider the "laminar motion" given by u ' = %' = 0, w Q f ~/cy', where k is 

a constant. Let the direction cosines of the axes of the ellipsoid (i.e., the 

.axes #, y, z) referred to the axes of x' t y' t z' be (l h m h n{), (l 2 , m 2) n 2 ) and 
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(h, wis, W) respectively. Then employing the usual formulae for the trans- 
formation of axes, we obtain 



"> 



j 



a = fcmiUi, b = fcm2?i2, C = tfm 3 %, 

f = \ K (m 2 nz + m 3 n 2 ), g = § ft (mtfij -f ?%%), h = J /e (mi% 4- m 2 fli), 

f = J k (m 2 % — m 3 ^ 2 ), v = h* ( m 3^i — mi %), f == i * ( m i% — wg?ii ), (40 > 
and equations (37) become 

(b 2 + (?) (oi = k (IPmtfiz — c% 3 % 2 ) 

(c? -\- a 2 ) (02 = K{a?m^n 1 — a 2 m 1 nz) y. (41) 

(a 2 +6 2 )&) 3 = fc(a 2 min2—b 2 m 2 ni) 

These may be regarded as the equations of motion of the ellipsoid and give 
its spins for any position of its axes. It may be seen from these that a 
possible type of motion is that in which any one of the three axes of the 
ellipsoid lies permanently along the axis of x\ i.e., perpendicular both to the 
direction of the velocity and to that of the velocity gradient in the undis- 
turbed motion, while the ellipsoid rotates about this axis with a variable spin.. 
A typical motion of this class is that for which 

h, mi, n i = 1, 0, ; h, m 2 , n 2 = 0, cos x> si* 1 X 5 
h, m> % = 0, —sin %, cos % ; ©1 = % ; co 2 = ©3 = 0. (42> 

Equations (41) then give 

(b 2 -f c 2 ) x = * Q> 2 cos2 X + ° 2 s * n2 X)> 
oi", omitting a constant of integration which is obviously immaterial, 

tan x = - tan -|L . (43) 

The equations readily admit of integration in the case of an ellipsoid of 
revolution. Let b = c, and introduce the three Euler angles, the angle 
between the axes of x and x\ <f> the angle between the planes of x'y' and x f x r 
-ty the angle between the planes of x'x and xy. Then 

• • * * 

co! = $ cos 0-h ty, W2 — sin ^r — <j> sin cos ty, ©3 = cos ijr+ $ sin 5 sin ^ 

and 

mi = sin cos 0, ii\ = sin sin <£ 

m 2 = — sin sin ^ + cos cos <£ cos i|r 

n 2 = cos $ sin\|r 4- cos # sin <£ cos ty y. (44) 

m 3 = —sin <j> cos i/r— cos cos <£ sin ^ 

m = cos <£ cos -\|r — cos # sin $ sin i|r ^ 
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Substituting in (41), we have after some reduction, 

m m 

<j[> COS + i/r = -L/CCOS0, (45) 

(a 2 -\-b 2 )6 = k {a? — b 2 ) sin 6 cos 6 sin $ cos <j>, (46) 

(a 2 + b 2 )<j) = k (a 2 cos 2 <£ -f 6 2 sin 2 <£). (47) 

.From (47) we have, if £ = when <£ = 0, 

tan ^ = 2 tan J^L, (48) 

and dividing (46) by (47) and integrating, 

tan2 G == tT7~2 aji ■ *a - 2jl > ( 49 ) 

7<r (a 2 cos 2 <j!> + 6 J sm j </> v 7 

where & is a constant of integration. 

It appears that the motion of the ellipsoid is periodic. Its axis of revolu- 
tion describes a cone about the perpendicular to the plane of the undisturbed 
motion, to which its inclination varies between the limits tdirr 1 {ajk) and 
tan" 1 (&/&). 

This statement requires modification in the degenerate cases of a disc and a 
thin rod respectively. Putting a = in equations (45), (46), and (47), we 
find that a disc can take any position in which its faces are entirely composed 
of stream-lines of the undisturbed motion of the fluid. It thus moves through 

o 

the fluid edge-on with its axis at any inclination to the fluid shear. If we 
put b = c = in the same equations we find what is otherwise obvious, that a 
thin rod can set itself in any direction in a plane which is a plane of constant 
velocity for the undisturbed motion of the fluid. In what follows, however, 
we shall regard rods and discs as spheroids of small but finite thickness. 

This investigation reveals no tendency on the part of the ellipsoid to set its 
axis in any particular direction with regard to the undisturbed motion of the 
fluid. The motions corresponding to all values of h are consistent with the 
boundary conditions and with the approximated equations of motion. 

| 5. The Dissipation of Energy. 

As the velocity of the fluid is finite at great distances from the ellipsoid, 
the total rate of dissipation of energy will be infinite. We will accordingly 
calculate the excess of the rate of dissipation of energy over that which obtains 
in the same motion of the fluid in the absence of the ellipsoidal particle. 
Difficulties which arise in connection with the boundary conditions, make it 
imperative that we should define precisely the two fluid motions which we 
are thus comparing. 

Describe a sphere S, of large but finite radius E, whose centre is at the 
centre of the particle. Let the velocity of the fluid at points on S be given 
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by equations (2), and let the ellipsoid move freely under the fluid pressures on 
its surface so that relations (38) are satisfied. We will calculate the excess 
of the rate of dissipation of energy over that which would occur if the ellipsoid 
were removed, the sphere S completely filled with fluid, and the velocity at 
all points on its surface the same as in the first case. 

The direct calculation of the rate of dissipation of energy presents very 
great difficulties. For our purpose, it will be sufficient to consider the rate 
at which work must be done over the boundary of the sphere S in order to 
maintain the motion. This, of course, differs from the rate of dissipation of 
energy by the rate of increase of the kinetic energy of the fluid. For the 
periodic motions investigated in the last section, however, the average rate of 
work of maintaining the motion must be the same as the average rate of 
dissipation of energy. 

As we shall ultimately suppose the radius of the sphere S to tend to 
infinity, we may adopt such a degree of approximation that negative powers 
of E are neglected in the final result. This implies that the stresses must be 
correct to order E~ 3 and the velocities to order E~ 2 . 

At great distances from the origin we have, neglecting terms of order r~' s , 

ft = -I l/r, X3 = %2 = %3 = 0. (50) 

Substituting in (18), (19), (20), and remembering that (38) are satisfied, and 

writing 

q> = Arf + Btf + G& + 2¥yz + 2Gzx + 2Bn!f, (51) 

we have, neglecting terms of order r~ 4 , 

to = u — 4#<J>/r 5 , v = vo—4y<£>/r b , w = iu —4z&/r b , (52) 

while (21) gives 

p=p — 8{i<&/r b . (53) 

Superpose a motion, finite at all points within the sphere S, and such as 
to cancel the second terms in the right-hand members of equations (52) at 
points on its surface. Such a motion is easily found by the ordinary spherical 
harmonic methods, and when combined with the motion represented by 

(52), we have, 



4x (E 5 - r 5 ) _ , 5 (E 2 - r 2 ) 3<I> •> 
4y(B»-r«) 5 (E 2 ~r 2 ) 3<I> ^ 

°0 E 5 r 5 ^+ E 5 dy f> 



(54) 



42(E 5 -r 5 ) , , 5(E 2 -r 2 )3<l> 

EV 5 ^ E 5 dz 

which give on the surface of S 

p = p — 50^/E 5 . (55) 
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These velocities will not accurately satisfy the boundary conditions on the 
surface of the ellipsoid, but proceeding by successive approximations it is 
easily seen that the additional terms will be of order 11" 7 on the surface of S, 
and may therefore be neglected. The stresses at r = R are calculated from 
(54) and (55), by the aid of (32). Typical results are 

" 5 ^ , 4x 2 ^ 2x c& 



xx = — po + 2/xa-f 10/>6 ^ 1T -$+ — <J> — 



l_fi 5 R 7 R 5 3: 



m 



y~ - Zfil + lU/x ^— 7 cp r6 ^ R6 & j* . (56) 

From these we have 

xr = (xxx-\-yxy-\-zxzjr)~R 

= _|, + | ( a , + h , + gs) + l(V {g*-^}, (57) 

and two similar equations. 

The rate at which work is done in maintaining the motion is given by 

^ (icq xr 4- Vo yr 4- u\) zr) dS, 

where the integration extends over the surface of S. On substitution from, 
(57), this integral is readily evaluated and it gives 

f ttE 3 . 2 / x(a 2 + b 2 + c 2 + 2f 2 +2g 2 -j-2h 2 ) 

+ ^7r/*(aA + bB + cC + 2fF + 2gG + 2hH). (58) 

The first term is the work necessary to maintain the motion in the absence 
of the particle ; the second represents the addition due to its presence. On 
making the radius of the sphere S tend to infinity, this result becomes exact. 
On substitution from (31) and (39) this gives for the additional rate of doing 
work owing to the presence of the particle 

<*W 16 f «o // a 2 +/3o // b 2 + 7 o ,/ c 2 

I, n , ... n — '— ' _*L_M. fTf /A *< hi m ■■ m m n u . 111 . 11 . 1 nr 1 1 !L— ■ — ■ mmm mm . ir n. nn 1 

dt 3 ^ L/3o / V / + 7o // «o // + VA // 

2f 2 2g 2 2h 2 "\ r 

+ « ' (W + <?y fa' (# + a*y 7oV + 6 2 )J # 
It is interesting to note that if due regard is not paid to the conditions on 
the bounding sphere, and if the calculation is made directly from equations 
(52), the result is only one-fifth of that given above. 

In the case of an ellipsoid of revolution (6 = 6) we have 70', 70" respectively 

equal to #/, A"> an d it is easy to show that ^o // + 2a f/ = 2& 2 a '. Equation 

1? 
'59) then takes the simplified form 

dW = u> f "<>"&' | b 2 + c 3 + 2PJ(g 8 + h 8 ) l , fim 

W r ** 126*00' £0" 26 2 « ' /3o'0 2 +& 2 )J' { } 
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If the fluid is moving in laminar motion, this gives on substitution from 



(40) and (44) 
rlW 



dt 



^TTfltt 2 



«0 



// 



It ' Ol« 



1 

+ ~ ; COS 2 9 + 

0* 



sin 4 # sin 2 2<£ 
2 



sin 2 $ 



(81 



The average values of cos 2 # and sin 4 6 sin 2 2<£ are calculated from (48) and 
(49) as follows : 

7/.2 

rZ0 



##& 



2ir{a 2 + b 2 ) 



'2tt $£ 

COS 2 -yj ^0 

" 



v/{(a 2 + P)(& 2 + /v 2 )} 



and 



KObh 



2a 2 b 2 



sin 4 # sin 2 26 -7- d6 = 7-5 — 
Jo ^ dcf) T {a 2 — 



J. 



a 2 + b 2 + 2P 



27r(a 3 + & 2 )Jo^" r # ~ r ~ (« 2 -^) 2 U/{(a 2 + & 2 )(& 2 + & 2 )} 

If D is the average rate of dissipation of energy, (61) then gives 

2a% 2 f a 2 + b 2 +2k 2 



9 I 



I) 



|- 7r/i« 2 



.(«» _ J2)2 l^/ { ( ft J + p) (£2 + / c 2) } 



• "» : I— 



¥ 



f *o" , 1 
L2J s «o'/3 " 2JW 
1 2 



/«p'(a 2 +& 2 ) 
2 



} 



+ 



^{(a 2 + k 2 ){b 2 + k 2 )}\bW f3 '(a 2 + b 2 )j l i3,'{a 2 + b 2 )l (62) 

We will show that this has no maxima or minima for finite, non-zero 
values of k, Differentiating with regard to k 2 , we obtain 

dD 27r^ 2 [k 2 j3o /f {2b 2 «o'~(a 2 + h^ .„„. 

d{k?) W^'^'^' {(a 2 + k 2 ){b 2 + k 2 )Y! 2 * K } 

From their definition in equations (10), we see that V> V, /V> A>" are 
all positive, while #o" —*$>" is of the same sign as a 2 — b 2 . Further 

(b 2 -X)d\ 



25V-(a 2 + 5 2 )/3 / = (a 2 -b 2 ) 



f (a 2 — 6 2 ) 



o (a 2 + \fl 2 (b 2 -tXf 
\d\ 



o (« 2 4-\) 5 / 2 (& 2 + A,f 

after an integration by parts. The latter integral is clearly positive, and 

hence the left-hand side is of the same sign as a 2 —b 2 . It follows that the 

right-hand side of (63) never vanishes for finite values of k, and is of opposite 

sign to a 2 —b 2 . Accordingly, the motion which gives minimum average 

dissipation of energy will correspond to k = oo for a prolate spheroid, and to 

k = for an oblate spheroid. The average rate, of dissipation of energy will, 

in all cases, lie between the values corresponding to k = and k = oo . 

We will write 

D = Wjjut 2 (64) 



Particles Immersed in a Viscous Fluid. 



175 



where v is a factor and v = 4:7rab 2 f3, and is therefore the volume of the 
particle. Equation (62) then gives, — f or k = oo , 

v = l/a&W, (65) 



and for & = 



z^ 



a$ 3 <x 

77"+"^ 7 + 



(66) 



In the case of spheroids the integrals in (9) and (10) can be evaluated in 
finite terms. "We obtain, 

(1) for prolate spheroids, b = c = a cos 

aW = icosec 4 0{(2-5cos 3 0) see 4 + 3©}, 
a 5 /3 ' = cosee 4 <9 {2 + sec 2 (9-30}, 
a 3 /3 ,/ = cosec 4 <9{(2 + cos 2 <9)©-3}, 



(67) 



where 



© = cosec log tan ( - -\- ~ ), 



(2) for oblate spheroids, b = c = a sec #, 

&<V = i : cosec 4 {30 -cos (5 -2 cos 2 0)}, 

&p ' = cosec 4 (9 {(2 + cos 2 0) sec (9 -3©}, 

jsfl," = cosec 4 (9 {(1 + 2 cos 2 (9) 0-3 cos 0} t (68) 

where © = #/sin 0. 

The maximum and minimum values of v for spheroids of different shapes 
are given in Tables I and II. They are given in terms of the " ellipticity " 
of the meridian section of « the particle, as this is generally more convenient 
for measurement than the eccentricity. The ellipticity e is defined to be the 
difference of the greatest and least diameters divided by the greatest, so that 
e = {a — b)ja for a prolate spheroid and e = (b — a)/b for an oblate spheroid. 



Table I. — Prolate Spheroids. 



Table II.— Oblate Spheroids. 



Ellipticity 
€ = « — b 

a 


Minimum 


Maximum 
v. 




O'l 

0*2 

0*3 

0*4 

0-5 

0-6 

0-7 

0-8 

0*9 

l'O 


2-5 

2*431 

2-361 

2-295 

2-232 

2-174 

2-120 

2*073 

2-035 

2-010 

2*000 


2*5 

2-540 
2-586 
2 -645 
2-719 
2-819 
2-958 
3-170 
3-548 
4-485 

00 



Ellipticity 
e = b~ a 


Minimum 
v. 


Maximum 
v. 


b ' 




0*1 

0-2 

0*3 

0-4 

0-5 

0-6 

0-7 

0-8 

0-9 

1-0 


2-5 

2*464 
2-426 
2-388 
2-348 
2*306 
2-262 
2-216 
2-168 
2-116 
2-061 


2-5 

2*582 
2-683 
2-818 
3*003 
3-267 

3 -670 

4 -354 
5-744 
9*960 

00 



vol. on. — A. 



N 
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It will be noted that the maximum value of v tends to infinity in each ease- 
as the ellipticity tends to unity. This does not, of course, mean that the 
dissipation of energy becomes infinite, since the volume of the particle at the 
same time tends to zero. 

| 6. The Viscosity of Suspensions. 

The results of the last section may be applied, as Einstein has applied the 
corresponding results for a sphere, to the determination of the viscosity of a 
fluid containing solid spheroidal particles in suspension. If the suspension is 
so dilute that the distances between neighbouring particles are great compared 
with their dimensions, we may regard each particle as the cause of a certain 
increase in the dissipation of energy. 

In the case of spherical particles the motion is steady and D is given 
accurately, without averaging, by the second term in (58). It is then easy to* 
show from (10), (25) and (39), that, when a = b = c } 

A/a = B/b = C/C = F/f = G/g = H/h = 5a 3 /8. (69) 

If v is the volume of the particle, so that v = 47m 3 / 3, we have, 

D = 5^(a 2 + b 2 + C 2 + 2t 2 +2g 2 + 2h 2 ). 

Hence, if v is the total volume of the particles in unit volume of the 
suspension, the total rate of dissipation of energy is 

2fi (1 + 2-5 V) (a 2 + b 2 + C 2 + W + 2g 2 + 2h 2 ). 

From the point of view of 'the dissipation of energy, we may, therefore,, 
regard the suspension as a homogeneous fluid whose coefficient of viscosity 

is jjl * where 

/** = /*(l + 2'5V). (70) 

This agrees with Einstein's well-known result. It does not require the 
particles to be equal in size, so long as they are all spherical in shape. It 
applies only to suspensions which are sufficiently dilute and in which there is 
no tendency to form aggregates. With these restrictions it is true for. any 
steady motion of the suspension, subject only to the condition that, if a, b, ... 
vary in space, they shall do so only on a scale which is large compared with 
the dimensions of the particles, and that this variation shall not promote a 
variation in the concentration of the particles. For example, in applying this 
formula to measurements of viscosity by capillary tube methods, we should 
have to be assured : (1) that the diameter of the tube is large compared with 
that of the particles ; (2) that there is no tendency for the particles to select 
particular paths, say, along the axis of the tube. 

Some further light may be thrown upon this last point. It is easily seen 
that, within the limitations of the present theory, there is no resultant force 
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oil the particle in a direction perpendicular to the axis of the tube. If the 
axis of the tube is taken as axis of z, the undisturbed motion is given by 
u = v = 0, w = A {a 2 — x 2 -*• 2/ 2 ) ? where a is the radius of the capillary tube 
and A is a known constant. Eeferred to parallel axes through the centre of 
a particle (h, Jc) t this gives u = v = ■ 0, w = & 2 — h 2 — & 2 — & 2 — £/ 2 — 2te — 2%. 
This may be regarded as the superposition of two motions : (1) u = v = 0, 
w = a 2 — h 2 — Jc 2 ■— x 2 — y 2 ; (2) it = # = 0, w = — 2te — 2%. Neither 
of these motions can give rise to any force on the particle perpendicular to 
the axis of z — the motion (1), for reasons of symmetry, and the motion 
(2), by the investigations of the present paper. It follows that the combined 
motion can give rise to no such force. 

In the case of laminar motion, we can extend Einstein's result to ellipsoidal 
particles. We have then from (64) for the average total dissipation of energy 
per unit volume of the suspension 

pfc 2 (l + vY) 
giving in place of (70) 

fi* = fi(l + vY). (71) 

The values of v for prolate and oblate spheroids of different ellipticity of 
meridian section are given in Tables I and II respectively. 

It will be noted that the minimum value of v is always less than the value 
(2*5) for spherical particles. In this connection it may be remarked that, 
although in some cases Einstein's formula agrees well with experiment, in 
others it agrees only if v has a value greater than 2*5.* Smoluchowskif 
proposes a formula of the type (71) for non-spherical particles, and, with 
regard to the factor v, he remarks: "For the sphere this value (v == 2*5) is 
evidently a minimum, since by virtue of its rolling motion it influences the 
shearing motion of the least volume of the fluid." This argument is not 
valid, for we have shown that for any spheroid there are possible motions for 
which the increase in the dissipation of energy is less than it is for a sphere 
of the same volume. 

§ 7. A Hypothesis on the Dissipation of Energy. 

It appears from the foregoing paragraphs that we are able only to specify 
limits for the numerical factor in Einstein's formula. It is obviously unde- 
sirable to leave a problem, which is physically quite determinate, in this 
indeterminate form. 

* For an account of the many experimental investigations of this problem, see 
Arrhenius, 'Biochem. Jour.,' vol. 11, p. 112 (1917). 

t " Theoretische Bemerkungen iiber die Viskositat der Kolloi'de," ' Kolloid Zeitsclk, 3 
vol. 18, p. 191 (1916). 
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We might assume that at any instant the axes of the particles are in 
random directions. On this hypothesis it is easy to calculate the rate of 
dissipation of energy. It is found that this gives a value for Einstein's 
factor which is greater than the maximum values shown in Tables I and II. 
The explanation is that random directions for the axes of the particles do 
not correspond to a permanent steady state. 

Again, we can determine v if we know the number of particles for which 
the k of (49) lies between k and Jc + dk, say f(k)dk. This, however, is just 
what our present theory fails to reveal: according to it, all values of k are 
possible, and there is no indication of their relative probability. It would 
appear that this failure is due to the limitations of the theory of the slow 
motion of a viscous fluid. 

It will be remembered that the only stable direction of permanent transla- 
tion of an ellipsoid in a perfect fluid is that of its least axis. If an attempt 
is made to solve the corresponding problem for a viscous fluid by Overbeck's 
method, it is found that no couple is exerted upon the ellipsoid by the fluid 
if its motion is one of pure translation in any direction whatever. A heavy 
ellipsoid falling through a viscous fluid would fall with different velocities for 
different directions of its axes with regard to the vertical, but apparently it 
would not tend to set its axes in any particular directions. This is due to 
the fact that in the case of the perfect fluid the resultant couples arise from 
those terms in the pressure which depend upon the square of the velocity, 
terms which are neglected in the ordinary viscous theory. 

It seems not improbable that, in our present problem, a more complete 
investigation would reveal the fact that the particles do tend to adopt special 
orientations with respect to the motion of the surrounding fluid. Our 
suggestion is that the particles will tend to adopt that motion which, of all the 
motions possible tmder the approximated equations, corresponds to the least 
dissipation of energy. If this be true, it reconciles the perfect fluid theory, 
the viscous theory, and experiment in the case of an ellipsoid falling under 
gravity through a fluid. 

Our assumption is reminiscent of the theorems of Helmholtz and Korteweg 
on the rate of dissipation of energy in the motion of a viscous fluid. A 
little reflection, however, will show that these theorems are not really 
applicable to our present problem. For one thing, except in the ease of 
spheroidal particles, no steady motion is possible. Further, these theorems 
are themselves deduced from the approximate equations of motion obtained 
by neglecting the squares and products of the velocities, and our investiga- 
tions have been based upon exact solutions of these same approximate 
equations. It is impossible that a general theorem, based upon a set of 
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differential equations, should give, with regard to a particular ease, informa- 
tion which is not contained in the corresponding exact solution of the 
equations. 

If this hypothesis be granted, we are able to assert the following 
propositions : — 

1. Prolate spheroidal particles immersed in a fluid in laminar motion will 
tend to set themselves with their axes perpendicular to the plane of the 
undisturbed motion of the fluid. They will rotate about their axes with 
constant angular velocity, and the whole motion will be steady. 

2. Oblate spheroidal particles immersed in a fluid in laminar motion will 
tend to set themselves with an equatorial diameter perpendicular to the plane 
of the undisturbed motion of the fluid. They will rotate about this diameter 
with a variable angular velocity, and the whole motion will be periodic but 
not steady. 

3. For motion in a capillary tube, the particles will tend to concentrate 
along the axis of the tube. 

4. The effective value of the numerical factor in Einstein's formula will be 
the minimum v, as shown in Tables I and II. This is always less than the 
2*5 appropriate to spherical particles. Thus, for particles of giv&n volume, 
the increase of the viscosity is greatest when the particles are spherical. 

An experimental investigation of some of these provisional conclusions 
would be valuable. 



